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In  this  paper  we  apply  the  theory  of  singularities  of  differentiable 
mappings  - specifically  the  unfolding  theorem  - to  study  the  effect  of 
imperfections  in  a system  subject  to  bifurcation.  In  a number  of  special 
cases  we  have  classified  (up  to  a suitable  equivalence)  all  the  possible 
perturbations  of  the  bifurcation  equations  by  a finite  number  of  imperfection 
parameters.  These  cases  include  both  bifurcation  from  a double  eigenvalue 
and  from  a simple  eigenvalue  degenerate  in  the  sense  of  Crandall-Rabinowitz . 
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SIGNIFICANCE  AND  EXPLANATION 


Bifurcation  theory  is  concerned  with  the  appearance  of  multiple  solutions  to 
the  equations  governing  the  equilibria  of  a physical  system  as  a parameter  is  varied. 
The  most  familiar  example  of  this  theory  concerns  the  buckling  of  a beam  as  the  com- 
pressive load  is  increased  (so-called  Euler  buckling)  - for  small  loads  the  unde- 
flected state  is  the  only  solution  of  the  governing  equations;  for  larger  loads 
solutions  with  deflection  also  exist,  and  are  in  fact  preferred  by  physical  systems, 
from  considerations  of  stability. 

It  is  widely  recognized  that  small  imperfections  in  a system  may  have  a profound 
influence  on  its  bifurcations.  For  instance,  in  the  idealized  treatment  of  the  beam 
problem  above,  the  theory  predicts  that  upwards  and  downwards  deflections  are  equally 
probable,  but  in  experiments  specimens  will  buckle  repeatedly  in  the  same  direction 
because  of  small  imperfections  in  the  experiment  which  destroy  the  symmetry  between 
up  and  down.  Also  the  experimentally  measured  buckling  loads  for  beams  are  dramati- 
cally less  than  predicted  by  the  idealized  theory.  Although  many  attempts  have  been 
made  to  explain  its  phenomenon,  no  universally  accepted  explanation  has  so  far  been 
produced . 

It  is  well  known  that  sensitivity  to  imperfections  is  greatly  increased  when  a 
system  loses  stability  in  situations  where  there  are  two  buckling  modes  present 
simultaneously  - i.e.  bifurcation  from  a double  eigenvalue. 

In  this  paper  we  use  the  theory  of  singularities  of  differentiable  mappings  to 
study  bifurcation  in  the  presence  of  imperfections.  In  a number  of  important  special 
cases,  including  some  bifurcations  from  a double  eigenvalue  and  some  bifurcations 
from  a simple  eigenvalue  degenerate  in  the  sense  of  Crandall-Rabinowitz , we  have 
classified  all  the  possible  perturbations  of  the  given  problem.  The  classification 
depends  on  a (small)  finite  number  of  imperfection  parameters. 

The  most  promising  applications  of  the  theory  concern  bifurcation  from  a double 
eigenvalue.  We  have  already  analyzed  bifurcations  near  a double  eigenvalue  of  the 
equations  governing  certain  model  chemical  reactions  with  these  techniques,  and  we 
anticipate  other  applications,  specifically  to  mode  jumping  in  buckled  plates  and 
to  the  selection  of  the  number  of  cells  in  the  Taylor  fluid  instability  problem  in 
a cylinder  of  finite  length. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  authors  of  this  report. 


equivalent  to  (1.4).  Finally,  the  theory 


A small  perturbation  of  one  sign  yields  a smooth 


of  the  attacks  that  have  been  directed  at  castrophe  theory  applications  of  catastrophe  theory  in  that  we  do  not  restrict  our  considerations 


which  can  occur  in  a problem.  Chow,  Hale,  and  Mallet-Paret  [ 6 ] have 


first  of  our  two  principal  theoretical  results  is  the  following 


LEMMA  2.7.  Suppose  G(x,X)  = xm  - Ax  and  H(x,A)  = any  other  problem  whose  Taylor  expansion  at  the  origin  agrees  with 


Finally  note  that  e(0)  = 1 since  H(x,0)  = xm.  diagrams  that  can  occur  in  the  universal  unfolding  F(x,A,a)  of 

The  next  proposition  gives  a more  theoretical  criterion  for  a given  problem.  In  this  discussion  we  assume  m * n.  As  mentioned 

when  two  bifurcation  problems  are  contact  equivalent.  A bifurca-  in  the  introduction*  there  are  three  ways  a diagram  can  fail  to  be 

tion  problem  G:  (£n  x r,o)  (Jgm,0)  is  called  k-determined  if  stable,  which  are  illustrated  in  Figure  2.1.  The  first  of  these, 


the  most  obviously  unstable  - a small  perturbation  will  in 


since  the  zero  set  of  F is  not  a manifold,  we  must  have  rank  dF 


but  not  a.  Thus  a diagram  with  bifurcation  can  occur  only  if 


germs  of  mappings.  Thus  for  example  the  bifurcation  problem 


1 


then  (3.9)  would  imply  (3.7).  Ke  now  find  such 


Section  4 Computations  with  one  degree  of  freedom  Proof:  Whether  or  not  (dll) (0)  is  zero  is  an  invariant  of 


When  m=2  and  index  (d  H) (0)  f 1 one  obtains  the  problem  y ♦ a • 

Such  problems  appear  in  the  chemical  engineering  literature  under  the  0ne  can  generalize  (III)  to  show  that  x3  - Akx  has 

name  of  isolas  121 ] • codimension  3k  * 1. 


thus  useful  to  be  able  to  determine  when  a 2-parameter  unfolding 


unfolding  of  G(x,A)  * x - (2A+l)sin  x.  This  is  the 


hysteresis  point  for  this  bifurcation  problem 


(x,X)  « x + a x*  ♦ 6x  - X is  the  universal 


(Ill)  and  (IV) . The  formalization  of  this  method  yields 


stituting  (4.10)  yields 


question  then  is  whether  the  bifurcation  problems  associated  to 
each  component  (in  the  algebraic  sense)  of  the  control  set  are 


all  different  for  8 near  0 so  we  have  proved  the  following: 


has  at  most  one  horizontal  and  two  vertical  tangents.  So  sequences 

ersal  unfolding  of  G.  The  following  relation  seems  to  hold  though  we  have 

_ ..  , . j like  LRRL  are  not  possible.  If  one  enumerates  all  the  possibilities 

codim G = No.  of  defining  conditions  -n  + modality.  r 

subject  to  the  above  constraints  one  finds  S3  distinct  stable 


in  Proposition  4.2  the  problem  G(x#A)  ■ x^h(x,A) 


geometrically  that  small  perturbations  of  this 


DEFINITION  5.2:  The  bifurcation  problem  (5.1)  is  nondegenerate 


also  note  that  the  Me Leod-Sattinger  results  [18]  follow  directly  from  Propo- 


Proof:  We  sketch  the  needed  calculations.  Writing  C = (g.h),  we 


and  p (z\  is  homogeneous  of  degree  two  In  x and  y.  Then 


• pi  it  into  a union  of  five  distinct  orbits,  described  by  the  structure  of  their 


by  (5. 14),  this  happens  only  for  finitely  many  w.  Figure  5.1  indicates  a domain 


quadratic  form,  is  part  of  the  standard  terminology  of  singularity  theory  [11  ]. 


solution  must  be  stable  for  some  range  of  the  parameter  and  therefore  has 


c* 

X 


Cfi 

O 

"O 


o 


c 


>» 

* 


C 3 


•-  00 


< 

lit 

w 4) 

£ & 


£ l 


l 

J9 


C _ 


"O 

£ 


3 

CO 

£ 


J:  ° 


o 

JC 


•5 

*«• 

o 

« 

c 

o 

« 

o 

CL 

*9 

a 


5 

5 


a £ 

o .2 


H 


8 ; 

3 I 

U O 

u Q. 

9 *-< 

!S  « 

a .2 

* • 

.2  *o 

« 


O 

a 


£ 

V 

3 

o 

ha 

a 

a 

o 


£ £ 


* 

o 

JC 

*3 

1 

o 


o 

la 

a 

o 

H 


t-  «• 

« ». 

O 5 


3 

9* 

*5 

o 

A 

v 

u 

n 

C 

O 

•9 


•9 

S 


c 

I 


the  finite  element  analogue  of  the  Euler 


be  written  symbolically  as  an  integration  by  parts,  when  permitted. 


6.21.  The  relations  (6.8)  are  satisfied.  second  column  contains  only  one  non-zero  element.  Expansion  in 


Our  strategy  in  evaluating  the  determinant  (6.9)  is  to  show  lemma  is  proved, 

the  last  row  contains  only  one  non-zero  element  and  that  the 


substituting  ♦Is)  - v Is)  - sin  s and  integrating  by  parts  universal  unfolding  for  the  idealised  bean,  problem. 


In  this  chapter  we  collect  a number  of  physical  problems  in 


states.  This  system  has  a bifurcation  diagram  of  the 


only  one  reactant  and  taking  place  in  a stirred  tank  in  which  the  (7.5) 


curved  plate  ( 23  , although  the  viewpoint  there  is  quite 


Postcn  and  Stewart's  forthcoming  book  ligl;  we  are  indebted  to  that  is,  only  if 


of  stability  for  the  branching  solutions 


higher  order  singularity  - specifically 


valent  function,  and  if  we  define  a new  variable  by 
[ (2X  -X)/3]  we  obtain  a new  function 


First,  let  us  consider  the  ODE  associated  with  space-independent 


Consider  a bifurcation  problem  in  variational  form. 
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